Large-Scale Virtual Screening for the Discovery of SARS-CoV-2 Papain-like Protease (PLpro) Non-covalent Inhibitors.

The rapid global spread of the SARS-CoV-2 virus facilitated the development of novel direct-acting antiviral agents (DAAs). The papain-like protease (PLpro) has been proposed as one of the major SARS-CoV-2 targets for DAAs due to its dual role in processing viral proteins and facilitating the host's immune suppression. This dual role makes identifying small molecules that can effectively neutralize SARS-CoV-2 PLpro activity a high-priority task. However, PLpro drug discovery faces a significant challenge due to the high mobility and induced-fit effects in the protease's active site. Herein, we virtually screened the ZINC20 database with Deep Docking (DD) to identify prospective noncovalent PLpro binders and combined ultra-large consensus docking with two pharmacophore (ph4)-filtering strategies. The analysis of active compounds revealed their somewhat-limited diversity, likely attributed to the induced-fit nature of PLpro's active site in the crystal structures, and therefore, the use of rigid docking protocols poses inherited limitations. The top hits were assessed against recombinant viral proteins and live viruses, demonstrating desirable inhibitory activities. The best compound VPC-300195 (IC50: 15 µM) ranks among the top noncovalent PLpro inhibitors discovered through in silico methodologies. In the search for novel SARS-CoV-2 PLpro-specific chemotypes, the identified inhibitors could serve as diverse templates for the development of effective noncovalent PLpro inhibitors.

Discovery of lead natural products for developing pan-SARS-CoV-2 therapeutics.

The COVID-19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), remains a global public health crisis. The reduced efficacy of therapeutic monoclonal antibodies against emerging SARS-CoV-2 variants of concern (VOCs), such as omicron BA.5 subvariants, has underlined the need to explore a novel spectrum of antivirals that are effective against existing and evolving SARS-CoV-2 VOCs. To address the need for novel therapeutic options, we applied cell-based high-content screening to a library of natural products (NPs) obtained from plants, fungi, bacteria, and marine sponges, which represent a considerable diversity of chemical scaffolds. The antiviral effect of 373 NPs was evaluated using the mNeonGreen (mNG) reporter SARS-CoV-2 virus in a lung epithelial cell line (Calu-3). The screening identified 26 NPs with half-maximal effective concentrations (EC50) below 50 µM against mNG-SARS-CoV-2; 16 of these had EC50 values below 10 µM and three NPs (holyrine A, alotaketal C, and bafilomycin D) had EC50 values in the nanomolar range. We demonstrated the pan-SARS-CoV-2 activity of these three lead antivirals against SARS-CoV-2 highly transmissible Omicron subvariants (BA.5, BA.2 and BA.1) and highly pathogenic Delta VOCs in human Calu-3 lung cells. Notably, holyrine A, alotaketal C, and bafilomycin D, are potent nanomolar inhibitors of SARS-CoV-2 Omicron subvariants BA.5 and BA.2. The pan-SARS-CoV-2 activity of alotaketal C [protein kinase C (PKC) activator] and bafilomycin D (V-ATPase inhibitor) suggest that these two NPs are acting as host-directed antivirals (HDAs). Future research should explore whether PKC regulation impacts human susceptibility to and the severity of SARS-CoV-2 infection, and it should confirm the important role of human V-ATPase in the VOC lifecycle. Interestingly, we observed a synergistic action of bafilomycin D and N-0385 (a highly potent inhibitor of human TMPRSS2 protease) against Omicron subvariant BA.2 in human Calu-3 lung cells, which suggests that these two highly potent HDAs are targeting two different mechanisms of SARS-CoV-2 entry. Overall, our study provides insight into the potential of NPs with highly diverse chemical structures as valuable inspirational starting points for developing pan-SARS-CoV-2 therapeutics and for unravelling potential host factors and pathways regulating SARS-CoV-2 VOC infection including emerging omicron BA.5 subvariants.

X-ray crystallographic characterization of the SARS-CoV-2 main protease polyprotein cleavage sites essential for viral processing and maturation.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the pathogen that causes COVID-19, produces polyproteins 1a and 1ab that contain, respectively, 11 or 16 non-structural proteins (nsp). Nsp5 is the main protease (Mpro) responsible for cleavage at eleven positions along these polyproteins, including at its own N- and C-terminal boundaries, representing essential processing events for viral assembly and maturation. Using C-terminally substituted Mpro chimeras, we have determined X-ray crystallographic structures of Mpro in complex with 10 of its 11 viral cleavage sites, bound at full occupancy intermolecularly in trans, within the active site of either the native enzyme and/or a catalytic mutant (C145A). Capture of both acyl-enzyme intermediate and product-like complex forms of a P2(Leu) substrate in the native active site provides direct comparative characterization of these mechanistic steps as well as further informs the basis for enhanced product release of Mpro's own unique C-terminal P2(Phe) cleavage site to prevent autoinhibition. We characterize the underlying noncovalent interactions governing binding and specificity for this diverse set of substrates, showing remarkable plasticity for subsites beyond the anchoring P1(Gln)-P2(Leu/Val/Phe), representing together a near complete analysis of a multiprocessing viral protease. Collectively, these crystallographic snapshots provide valuable mechanistic and structural insights for antiviral therapeutic development.

Ligand Binding Prediction Using Protein Structure Graphs and Residual Graph Attention Networks.

Computational prediction of ligand-target interactions is a crucial part of modern drug discovery as it helps to bypass high costs and labor demands of in vitro and in vivo screening. As the wealth of bioactivity data accumulates, it provides opportunities for the development of deep learning (DL) models with increasing predictive powers. Conventionally, such models were either limited to the use of very simplified representations of proteins or ineffective voxelization of their 3D structures. Herein, we present the development of the PSG-BAR (Protein Structure Graph-Binding Affinity Regression) approach that utilizes 3D structural information of the proteins along with 2D graph representations of ligands. The method also introduces attention scores to selectively weight protein regions that are most important for ligand binding. Results: The developed approach demonstrates the state-of-the-art performance on several binding affinity benchmarking datasets. The attention-based pooling of protein graphs enables identification of surface residues as critical residues for protein-ligand binding. Finally, we validate our model predictions against an experimental assay on a viral main protease (Mpro)-the hallmark target of SARS-CoV-2 coronavirus.

Automated discovery of noncovalent inhibitors of SARS-CoV-2 main protease by consensus Deep Docking of 40 billion small molecules† † Electronic supplementary information (ESI) available. See DOI: 10.1039/d1sc05579h

Recent explosive growth of ‘make-on-demand’ chemical libraries brought unprecedented opportunities but also significant challenges to the field of computer-aided drug discovery. To address this expansion of the accessible chemical universe, molecular docking needs to accurately rank billions of chemical structures, calling for the development of automated hit-selecting protocols to minimize human intervention and error. Herein, we report the development of an artificial intelligence-driven virtual screening pipeline that utilizes Deep Docking with Autodock GPU, Glide SP, FRED, ICM and QuickVina2 programs to screen 40 billion molecules against SARS-CoV-2 main protease (Mpro). This campaign returned a significant number of experimentally confirmed inhibitors of Mpro enzyme, and also enabled to benchmark the performance of twenty-eight hit-selecting strategies of various degrees of stringency and automation. These findings provide new starting scaffolds for hit-to-lead optimization campaigns against Mpro and encourage the development of fully automated end-to-end drug discovery protocols integrating machine learning and human expertise. Deep learning-accelerated docking coupled with computational hit selection strategies enable the identification of inhibitors for the SARS-CoV-2 main protease from a chemical library of 40 billion small molecules.

Automated discovery of noncovalent inhibitors of SARS-CoV-2 main protease by consensus Deep Docking of 40 billion small molecules.

Recent explosive growth of 'make-on-demand' chemical libraries brought unprecedented opportunities but also significant challenges to the field of computer-aided drug discovery. To address this expansion of the accessible chemical universe, molecular docking needs to accurately rank billions of chemical structures, calling for the development of automated hit-selecting protocols to minimize human intervention and error. Herein, we report the development of an artificial intelligence-driven virtual screening pipeline that utilizes Deep Docking with Autodock GPU, Glide SP, FRED, ICM and QuickVina2 programs to screen 40 billion molecules against SARS-CoV-2 main protease (Mpro). This campaign returned a significant number of experimentally confirmed inhibitors of Mpro enzyme, and also enabled to benchmark the performance of twenty-eight hit-selecting strategies of various degrees of stringency and automation. These findings provide new starting scaffolds for hit-to-lead optimization campaigns against Mpro and encourage the development of fully automated end-to-end drug discovery protocols integrating machine learning and human expertise.

COVID-19 Knowledge Extractor (COKE): A Curated Repository of Drug-Target Associations Extracted from the CORD-19 Corpus of Scientific Publications on COVID-19.

The COVID-19 pandemic has catalyzed a widespread effort to identify drug candidates and biological targets of relevance to SARS-COV-2 infection, which resulted in large numbers of publications on this subject. We have built the COVID-19 Knowledge Extractor (COKE), a web application to extract, curate, and annotate essential drug-target relationships from the research literature on COVID-19. SciBiteAI ontological tagging of the COVID Open Research Data set (CORD-19), a repository of COVID-19 scientific publications, was employed to identify drug-target relationships. Entity identifiers were resolved through lookup routines using UniProt and DrugBank. A custom algorithm was used to identify co-occurrences of the target protein and drug terms, and confidence scores were calculated for each entity pair. COKE processing of the current CORD-19 database identified about 3000 drug-protein pairs, including 29 unique proteins and 500 investigational, experimental, and approved drugs. Some of these drugs are presently undergoing clinical trials for COVID-19. The COKE repository and web application can serve as a useful resource for drug repurposing against SARS-CoV-2. COKE is freely available at https://coke.mml.unc.edu/, and the code is available at https://github.com/DnlRKorn/CoKE.

Deep Modeling of Regulating Effects of Small Molecules on Longevity-Associated Genes.

Aging is considered an inevitable process that causes deleterious effects in the functioning and appearance of cells, tissues, and organs. Recent emergence of large-scale gene expression datasets and significant advances in machine learning techniques have enabled drug repurposing efforts in promoting longevity. In this work, we further developed our previous approach-DeepCOP, a quantitative chemogenomic model that predicts gene regulating effects, and extended its application across multiple cell lines presented in LINCS to predict aging gene regulating effects induced by small molecules. As a result, a quantitative chemogenomic Deep Model was trained using gene ontology labels, molecular fingerprints, and cell line descriptors to predict gene expression responses to chemical perturbations. Other state-of-the-art machine learning approaches were also evaluated as benchmarks. Among those, the deep neural network (DNN) classifier has top-ranked known drugs with beneficial effects on aging genes, and some of these drugs were previously shown to promote longevity, illustrating the potential utility of this methodology. These results further demonstrate the capability of "hybrid" chemogenomic models, incorporating quantitative descriptors from biomarkers to capture cell specific drug-gene interactions. Such models can therefore be used for discovering drugs with desired gene regulatory effects associated with longevity.

Comprehensive Consensus Analysis of SARS-CoV-2 Drug Repurposing Campaigns.

The current COVID-19 pandemic has elicited extensive repurposing efforts (both small and large scale) to rapidly identify COVID-19 treatments among approved drugs. Herein, we provide a literature review of large-scale SARS-CoV-2 antiviral drug repurposing efforts and highlight a marked lack of consistent potency reporting. This variability indicates the importance of standardizing best practices-including the use of relevant cell lines, viral isolates, and validated screening protocols. We further surveyed available biochemical and virtual screening studies against SARS-CoV-2 targets (Spike, ACE2, RdRp, PLpro, and Mpro) and discuss repurposing candidates exhibiting consistent activity across diverse, triaging assays and predictive models. Moreover, we examine repurposed drugs and their efficacy against COVID-19 and the outcomes of representative repurposed drugs in clinical trials. Finally, we propose a drug repurposing pipeline to encourage the implementation of standard methods to fast-track the discovery of candidates and to ensure reproducible results.

A critical overview of computational approaches employed for COVID-19 drug discovery.

COVID-19 has resulted in huge numbers of infections and deaths worldwide and brought the most severe disruptions to societies and economies since the Great Depression. Massive experimental and computational research effort to understand and characterize the disease and rapidly develop diagnostics, vaccines, and drugs has emerged in response to this devastating pandemic and more than 130 000 COVID-19-related research papers have been published in peer-reviewed journals or deposited in preprint servers. Much of the research effort has focused on the discovery of novel drug candidates or repurposing of existing drugs against COVID-19, and many such projects have been either exclusively computational or computer-aided experimental studies. Herein, we provide an expert overview of the key computational methods and their applications for the discovery of COVID-19 small-molecule therapeutics that have been reported in the research literature. We further outline that, after the first year the COVID-19 pandemic, it appears that drug repurposing has not produced rapid and global solutions. However, several known drugs have been used in the clinic to cure COVID-19 patients, and a few repurposed drugs continue to be considered in clinical trials, along with several novel clinical candidates. We posit that truly impactful computational tools must deliver actionable, experimentally testable hypotheses enabling the discovery of novel drugs and drug combinations, and that open science and rapid sharing of research results are critical to accelerate the development of novel, much needed therapeutics for COVID-19.

COVID-19 Knowledge Extractor (COKE): A Tool and a Web Portal to Extract Drug - Target Protein Associations from the CORD-19 Corpus of Scientific Publications on COVID-19.

Objective: The COVID-19 pandemic has catalyzed a widespread effort to identify drug candidates and biological targets of relevance to SARS-COV-2 infection, which resulted in large numbers of publications on this subject. We have built the COVID-19 Knowledge Extractor (COKE), a web application to extract, curate, and annotate essential drug-target relationships from the research literature on COVID-19 to assist drug repurposing efforts. Materials and Methods: SciBiteAI ontological tagging of the COVID Open Research Dataset (CORD-19), a repository of COVID-19 scientific publications, was employed to identify drug-target relationships. Entity identifiers were resolved through lookup routines using UniProt and DrugBank. A custom algorithm was used to identify co-occurrences of protein and drug terms, and confidence scores were calculated for each entity pair. Results: COKE processing of the current CORD-19 database identified about 3,000 drug-protein pairs, including 29 unique proteins and 500 investigational, experimental, and approved drugs. Some of these drugs are presently undergoing clinical trials for COVID-19. Discussion: The rapidly evolving situation concerning the COVID-19 pandemic has resulted in a dramatic growth of publications on this subject in a short period. These circumstances call for methods that can condense the literature into the key concepts and relationships necessary for insights into SARS-CoV-2 drug repurposing. Conclusion: The COKE repository and web application deliver key drug - target protein relationships to researchers studying SARS-CoV-2. COKE portal may provide comprehensive and critical information on studies concerning drug repurposing against COVID-19. COKE is freely available at https://coke.mml.unc.edu/ and the code is available at https://github.com/DnlRKorn/CoKE.

Crystallographic structure of wild-type SARS-CoV-2 main protease acyl-enzyme intermediate with physiological C-terminal autoprocessing site.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the pathogen that causes the disease COVID-19, produces replicase polyproteins 1a and 1ab that contain, respectively, 11 or 16 nonstructural proteins (nsp). Nsp5 is the main protease (Mpro) responsible for cleavage at eleven positions along these polyproteins, including at its own N- and C-terminal boundaries, representing essential processing events for subsequent viral assembly and maturation. We have determined X-ray crystallographic structures of this cysteine protease in its wild-type free active site state at 1.8 Å resolution, in its acyl-enzyme intermediate state with the native C-terminal autocleavage sequence at 1.95 Å resolution and in its product bound state at 2.0 Å resolution by employing an active site mutation (C145A). We characterize the stereochemical features of the acyl-enzyme intermediate including critical hydrogen bonding distances underlying catalysis in the Cys/His dyad and oxyanion hole. We also identify a highly ordered water molecule in a position compatible for a role as the deacylating nucleophile in the catalytic mechanism and characterize the binding groove conformational changes and dimerization interface that occur upon formation of the acyl-enzyme. Collectively, these crystallographic snapshots provide valuable mechanistic and structural insights for future antiviral therapeutic development including revised molecular docking strategies based on Mpro inhibition.

Rapid Identification of Potential Inhibitors of SARS-CoV-2 Main Protease by Deep Docking of 1.3†Billion Compounds.

The recently emerged 2019 Novel Coronavirus (SARS-CoV-2) and associated COVID-19 disease cause serious or even fatal respiratory tract infection and yet no approved therapeutics or effective treatment is currently available to effectively combat the outbreak. This urgent situation is pressing the world to respond with the development of novel vaccine or a small molecule therapeutics for SARS-CoV-2. Along these efforts, the structure of SARS-CoV-2 main protease (Mpro) has been rapidly resolved and made publicly available to facilitate global efforts to develop novel drug candidates. Recently, our group has developed a novel deep learning platform - Deep Docking (DD) which provides fast prediction of docking scores of Glide (or any other docking program) and, hence, enables structure-based virtual screening of billions of purchasable molecules in a short time. In the current study we applied DD to all 1.3 billion compounds from ZINC15 library to identify top 1,000 potential ligands for SARS-CoV-2 Mpro protein. The compounds are made publicly available for further characterization and development by scientific community.